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Abstract 

We analyze the renormalization of wave functionals and energy eigenvalues in field 
theory. A discussion of the structure of the renormalization group equation for a general 
Hamiltonian system is also given. 



1. Introduction 

The Schrodinger representation of a field theory, using wave functions which are func- 
tionals of the field variables, has been known for a long time [1]. Nevertheless, it has rarely- 
been used since other representations are far more convenient for most field theoretic cal- 
culations. However, there are many problems, especially those involving changes in ground 
state properties, for which our intuitive understanding is better served in a Hamiltonian 
framework and using wave functionals. Analysis of Quantum Chromodynamics in light- 
cone framework [2] , saturation problems in deep inelastic scattering [3] and many problems 
in condensed matter physics may belong to this category. Wave functionals may also be 
useful in a variational approach. The most obvious difficulty in using wave functionals in 
field theory is in understanding their renormalization properties, especially in the modern 
or Wilsonian interpretation, where we consider how various quantities must change as high 
momentum modes are progressively integrated out [4] . Actually, some time ago, motivated 
by application to the Casimir effect, Symanzik analyzed a closely related question [5]. He 
considered the Euclidean partition function in a space with a boundary and used com- 
binations of Dirichlet and Neumann-type Green's functions to ensure proper behaviour 
of field variables on the boundary. Renormalization can be carried out by introducing 
counterterms which are monomials of the field and their derivatives on the boundary, in 
addition to the usual counterterms in the action. Although the results can be reinterpreted 
using wave functions [5,6], the analysis is not in the context of partially integrating out 
high momentum modes. In this paper, we study the change of wave functionals as high 
momentum modes are integrated out and obtain the corresponding renormalization group 
(RG) equation. This may also be the first step in alleviating some of the difficulties with 
variational wave functions in field theory [7] . 
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In the next section we shall consider the theory of a real scalar field since this suffices 
to illustrate most of the general ideas. In section 3, we show explicitly the computation 
of the renormalization constants related to the surface counterterms. The RG-equations 
for wave functionals and energy eigenvalues are obtained in section 4. Although these 
equations are derived for our example of a scalar field, the structure is quite general and 
will hold in any field theory. The last section gives a discussion of the RG-equation for a 
general Hamiltonian system. 

2. The reduced wave functional 

Consider the theory of a real scalar field (f){x, t) with Z2 symmetry with an action of 
the form 

S{A,(P,t,tQ) = I dx^ I d^x Z3(A) [l{d4>f + l{m^ - 6m^)4>^] + Z^{A)X<f>^ (1) 

A denotes the upper cutoff on the momenta. The wave functions of the theory are func- 
tionals of the field at the fixed time t. We decompose the original field into high and low 
momentum modes as 0(a;, t) — x(x, t) + ip{x^ t), where x denotes modes of high momenta 
(heavy modes), i.e. with momenta constrained to a shell /j, < k < A, and (p denotes modes 
of low momenta (light modes) with k < /j,. 

Integration over the heavy modes generates an effective field theory with a cutoff /i. Of 
course, this integration entails also the redefinitions of parameters so that the low energy 
correlation functions are preserved, i.e., they are invariants of the renormalization group 
(RG) flow. In terms of the wave function, this must produce a low energy wave function 
'^redi'p) starting from '^{(j)) = ^[x,<p]- As we shall see below in explicit calculations, this 
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reduced wave function has the form 



'^red{l^,V,t) = U[n,(p] / [d(p] e*' 



(2) 



where U is given by 



U[iJ,,ip] =expi / d^x [Z4(iJ,)(fidtip =e 



(3) 



U compensates for possible cutoff dependent surface terms ( on the time-shce at t) gener- 
ated by the integrations over x-modes. (p) is given by Eq.(l) with a cutoff fj, instead 
of A. For the wave function, we need new renormahzation constants Z4 and S. The inte- 
gration in Eq.(2) is over aU values of the fields between to and t, and over the fields at to, 
as usual. 

Now, integration over the heavy modes, at least diagrammatically, means that we 
allow x-fields only on internal lines or propagators. Such propagators represent vacuum 
fluctuations of the x-field. Hence we can obtain the reduced wave function by treating the 
X-modes as a kind of environment to which the v'-modes are coupled, the distribution of 
the heavy modes in the environment being characterized by a vacuum wave function of 
the x-niodes, say ^PoM- Thus we start with \E'(A) which has a representation as in Eq.(2), 
i.e.. 





Writing = x + V', the action splits as 



(5) 



We define a function 



F{x,v,t)= [dx] expi[S{x,t,to) + Sint{x,<P,t,to)] *(x,^,to) (6) 
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The initial wave function can be assumed to factorize into a x-dependent term and a 
(/^-dependent term, at least as to ~^ —oo; this is essentially the Born-Oppenheimer approx- 
imation. ^o[x] = U[x,<p = 0] F{xi'P = 0) then defines the "vacuum wave function" for 
the x-modes. The integration over x-niodes can then be represented on the wave function 
by *(0,A) ^ *r-ed(¥',/i) with 

*red(<^,At) = J [dx] *o(x,A,t) *(x,<^,A,t) = < nx\x,V> (7) 

This formula defines the procedure of integrating out the high momentum modes and, 
along with (2), will lead to the computation of the renormalization constants Z4 and S. 

Consider now the time-evolution of ^red as given by Eq.(2). The quantity U[(f] 
involves variables on the time-slice at t and can be interpreted in the operator language 
as a unitary transformation. Keeping this in mind, the variation of the time-argument in 
Eq.(2) gives 

where 

Hred = U H{^) (9) 

and H{ip) is the Hamiltonian operator given by the action S{iJ,,(p,t,tQ); H(ip) includes 
some of the effects of integrating out the x's, viz. those associated with spacetime volume 
integrals. By rewriting '^red using Eq.(7), we see that 

idt'ifredif,!^) =< - Hx\x,<P> + < ^^xl^(x,^)lx,^ > (10) 

If the x-inodes have momenta very high compared to the <^'s of interest, they cannot be 
excited by the (p^s for kinematical reasons. We can expand \x, > in terms of a complete 
set of x-states and we then see that the above matrix elements vanish except for the case 
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Ix, ^ >= *o(x)l¥' >) X"n^odes have only vacuum fluctuations. Eq.(lO) thus gives -ffred 



as 



Hredi<fi,l^) =< ^x\Hix,<f)\^x > -Eq^ (11) 

If x's are thought of as environmental degrees of freedom, we expect the effective Hamilto- 
nian for (/p's to be H{x, averaged over the distribution of x's, exactly as in Eq.(ll). We 
now see that such a Hamiltonian does indeed generate the time-evolution of the reduced 
wave function as defined in Eq.(2). 

3. Explicit calculation of the reduced wave function 

We now proceed to calculate explicitly the reduced wave function as defined by Eq. (7) . 
There will be the usual corrections to the action which depend on the cutoff and which are 
spacetime volume integrals. Here we are primarily interested in the surface terms which 
contribute to U[ip\ and hence to Z4, 5. We shall calculate these to first order in A. 

The (^-modes play the role of background fields as far as the calculation of F(x, is 
concerned. From its definition, we can write down an equation for F(x, ^) as 



idtF 



(12) 

For calculations to first order in A, we can simplify this as 



I ld'x[-^+u;'x' + Vx']F (13) 



where a;^x^ = x(~V^ + m^)x and V = 12X(p^. V varies slowly in time compared to x- 
Thus an adiabatic approximation is suitable and we can use a WKB-type ansatz 

F = exp(z(5(°) + + ...)) (14) 
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with dtS^^^ C'(l), S^S^-^^ ^ 0{1) , etc. To the order we are interested in, we get 



F exp 



'to 

The reduced wave function has the form 



(15) 



*red = J [d<p] eMiS{<p)) J F* (x, <p = 0)F{x, <p)dx (16) 

By expanding F{x, <f) in a series in V we obtain a hnearly divergent surface term 

V 3A f 

and a logarithmicaUy divergent surface term 

1<J X^X >=i^ log(AV/.2) J ^dt^d'x, {17b) 

where < X >= J dx F*{xj 'P = 0)XF(x, (p = 0). CoUecting terms we get to first order in 
A 

iG{fi) = iG{A) +^{A-^)J d^x ^2 + il^ log(AV/i') J d^x pdtp + ... (18) 



For the A-dependence to cancel out we need 

^4(A) +01og(AV^^) =ZM 
3A 

i6{A) +^(A-//) ^i6{i^) 

or in other words, 

3A 

Z^{A) =C4 -^logA^ +... 



(19) 



(20) 



The constants C4, cq depend on the normahzation conditions chosen; a minimal choice is 
to set them to zero. (Symanzik used C4 = 1.) 
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One can give general power counting arguments, along the lines of ref. [4] , to see that 
the surface counterterms we have introduced suffice to absorb all the A-dependence. 



4. RG Equation for * 

The wave function (3) is independent of A, as seen from the explicit computations, or 
equivalently, ^ in Eq.(2) is independent of fx when further integrations are carried out. This 
leads to the renormalization group equation. In what follows, we shall be using reduced 
wave functions and Hamiltonians; we shall leave this qualification as understood, omitting 
the subscript red. Also we ignore the mass terms for simplicity; they can be restored in 
a straightforward fashion if desired. Note that Z^dtip is the canonical momentum n. We 
may thus write 

G = iZ^Z-^ J d^x (pn (21) 

Consider an infinitesimal change in the cutoff A in the action, say A — > A(l — e). After 
Fourier transformation and rescaling of the corresponding momenta, we find 



dS 



j d'^x ix-Vip + ip + ^ifi)-^S-/3~S + J Hdt (22) 



cZlogA 

The origin of the last term (J Hdt) is due to the fact that the integration over time is only 
upto t. In Eq.(22), we have introduced both the function and the anomalous dimension 
7, defined by 

1 

^=-2 



^logZa 



P = - 



9 log A 



(23) 



A(o) 



91og A_ 

where A^°) is the bare coupling constant, defined by Z^X = Z^X^^\ 

Let us now consider the surface term Z/^Z^^ / V'^- It is easy to see, again by Fourier 
transformation followed by rescaling of momenta, that J (f tt is independent of A; the 



J [dip] J i{x-V(p + (f + 'y(p)^ 



contribution to the variation of the wave function from G{A) is thus solely due to ^4^3 ^. 
In using Eq.(22) we need the partial integration formula 

exp{iS) — J d^x {x ■ Wip + (/? + ^(p)-^ ( j [dip] exp(i5')) 

(24) 

On the right hand side, the field variables and functional derivative refer to the values 
on the boundary, i.e., at time t. Using this equation and the A- independence of we 
immediately get the renormalization group equation 



dlogA 
or equivalently 



d^x (x • + (1 + 7 + a)ip)] TT - (3— + iH t 



^' = 0. (25a) 



d^x (x • + (1 + 7 + a)(p)TT - (3— - t 



d .d_ 

m 



* = 



(256) 



where 



cr = — 



9(^4^3-^ 
9 log A 



A(0) 



3A 
4^ 



+ 



(26) 



is the new anomalous dimension arising due to the logarithmically divergent surface term. 
This result agrees with Symanzik [4]. 

Strictly speaking, our discussion so far applies to the vacuum wave function. This is 
because we have neglected the nontrivial (/?- dependence and associated properties of the 
initial wave function in the functional integral formulae. When we discuss excited states, 
we must have a nonzero number of "particles" in the initial state or they have to be created 
at some time prior to t which corresponds to the insertion of field variables in the functional 
integral. Either way, the partial integration formulae are more subtle. Notice that it is 
U~^'^ which is evolved by e*"^. Thus Eq.(2) should more accurately read as 



*(t) = U{t) j [dip] e*^(*'*°) (C/-i(to)*(to)) 



(27) 



The cutofF-independence of the wave functions then leads to the equation 



(Vm) - U{t) / e'^ U-\to)iVmo) = 



(28) 



where we have defined 



/d 
d^x [a; • V<^ + (1 + 7 + a)(p\7r - (3— + iHt 



(29) 



The question is whether we can set V'^ to zero. Let us say that are conserved quantities 
in the theory which can be used to label the states, i.e., Qa^q = Qa^q- Applying V on 
this equation, it is easy to see that V'^ = is inconsistent unless \D, Qa\ = 0. For the 
spatial momentum Pj, which is conserved in the present case, we see that [D, Pj] = Pj and 
we cannot set to zero. From the conservation property, we have 



Pi 



dpi 



exp(-^m)p,^^^ = U{t) J e'^ U{tor'p 



1 d^pito) 



dpi 



(30) 



Combining this with Eq.(28), we have 



9* fa. 5* 



(31) 



It is consistent to set {V + c P^)^ to zero if we choose c = 1. More generally, the 



RG-equation reads 



=0 



(32a) 
(326) 



Specifically for the case of N particles of momenta , N , we have 



* = (33) 
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This gives the RG-flow of the wave function. For the special case of energy eigenstates, 
viz., for 

^ = e-''^^ i;[(fi], (34) 
Eq.(33) also gives the flow of energy eigenvalue oj as 

a; -a; =0 (35) 

The solution to this equation is given, for example for a one-particle state, in terms of an 
arbitrary function f{s) as 

a;= V^/(« + log(p//i)) (36) 

where du = A simple application of this result is to the one-electron state in the BCS 
theory of superconductivity for which we have (3 — —aX^ for the four-fermion coupling, 
with a = rnkp/'^T^'^ [8]. ( Positive a corresponds to asymptotic freedom, m is the electron 
mass and kp is the Fermi momentum.) In this case, / = log(p//x) -|- 1/aA, and if one 
assumes a mass gap, /(s) ^ for small s, so as to cancel the p-dependence, eq. (36) 
then leads to the well-known relation lu = /j, exp 

Our derivation of the RG-equations has been in the context of the scalar field theory; 
nevertheless, it is clear that the general structure of Eqs.(25), (33) and (35) will be obtained 
in any field theory. 

5. RG-equation in Classical and Quantum Dynamics 

The basic idea of the renormalization group, viz. progressively integrating out the 
high momentum modes is very general. It is of some interest therefore to bring out the 
structure of the RG-equations already in classical dynamics. Consider a Hamiltonian 
system with phase space variables (^* ,%*), where x's represent heavy modes or variables 
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weakly coupled to the "light modes" One can integrate out the heavy modes to obtain 
an effective Hamiltonian for the variables (A specific, purely classical, context in which 
this can be concretely realized is in statistical mechanics where a Liouville distribution of 
X's can be given and they can be integrated out.) The change of the effective Hamiltonian 
as we lower the cutoff A is what we are interested in. We write 

H[A{1 - e)] = H{A) - ej^ ^ H + eV (37) 

The equations of motion with this Hamiltonian are 

where ^ujijd^^ A d^^ = uj is the symplectic two-form. In quantum field theory, we seek 
to preserve the correlation functions as we lower the cutoff. In the present case, what is 
preserved is the set of solutions to the equations of motion. Thus one should be able to 
absorb the change due to F as a change of variables ^* = + ^'^^{0^ ^ change of 

parameters — > g'°' = g"' + and a change of scale for the time variable t ^ t' = 
t{l — ehn)- (hn should be 1 for relativistic theories with identical scaling of space and time 
variables. It can be different from 1 for nonrelativistic theories where the time variable is 
scaled differently from the spatial coordinates.) In other words, V should be such that Eq. 
(38) can be written as 



■(l + eM = (^-'(e',^')r 



dH 



(39) 



dt 

Expanding to first order in e and using Eq. (38), we get 

n + - - /'"fl + h^PM = (40) 

where Pj = (w~^)'^'^ for f — H,V. Pj is the Hamiltonian vector field corresponding 
to the phase space function /. Eq. (40) is the basic expression of the invariance under 
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RG-flow. For many cases of interest, such as our field theory example, are taken to 

be linear in ^\ i.e., = -D^^'^ = {D$,y where D^. are constants. We also have constant 
ujij with 

-u.^Dliu-')'^^ = D] (41) 

In this case, we can multiply Eq. (40) by Uji and simplify to obtain 

d 



Thus an RG equation for H would read 



dg°' d^'^ 



= (42) 



d(logA) ("i^'w-^'W-^^"" '''' 
The term involving {D^Y can be written as the Poisson bracket {D^H} where D is the 
canonical generator of the transformation. For example, D — J d?x tt{x • V + 1 + 7 + a)(p 
for the scalar field theory of Eq. (1). Eq. (43) can thus be written as 

+ {D, H} +P°'^ + hHH = (44) 



d(logA) ^ ' dg' 
At this stage, one can also go to the quantum version in the standard way by replacing 
Poisson brackets by commutators. The resulting equation can also be derived from the 
compatibility of Schrodinger equation and the RG-equation (25). 
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